Spin-boson dynamics beyond conventional perturbation theories 
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A novel approximation scheme is proposed to describe the dynamics of the spin-boson problem. 
Being nonperturbative in the coupling strength nor in the tunneling frequency, it gives reliable 
results over a wide regime of temperatures and coupling strength to the thermal environment for 
a large class of bath spectral densities. We use a path-integral approach and start from the exact 
solution for the two-level system population difference in the form of a generalized master equation 
(GME). Then, we approximate inter-blip and blip-sojourns interactions up to linear order, while 
retaining all intra-blip correlations to find the kernels entering the GME in analytical form. Our 
approximation scheme, which we call Weakly-Interacting Blip Approximation (WIBA), fully agrees 
with conventional perturbative approximations in the tunneling matrix element (Non-Interacting 
Blip Approximation) or in the system-bath coupling strength. 

PACS numbers: 05.30.-d, 03.65.Yz, 05.40.-a 



I. INTRODUCTION 

Many physical and chemical two-level systems (TLSs) 
suffer the influence of external environments which cause 
decoherence effect*^ 3 -. Electron and proton transfer re- 
action in condensed phases 4 ^, tunneling phenomena in 
condensed matter physics^, or two-level atoms in an 
optical cavity 8 * are some very well known situations sus- 
ceptible to lack of coherence. Dissipation processes are 
to be understood first, in order to provide, if required, 
further decoherence control schemes. An important foot- 
step to allow experimental investigation of decoherence 
mechanisms was the realization of micrometer-sized ob- 
jects like the radio-frequency superconducting quantum 
interference device (rf-SQUID)£. Recently, due to their 
scalability and to the ease to experimentally control their 
parameters, other superconducting devices as the flux- 
qubil£ or the Cooper-pair bo x 10 : 11 have become attrac- 
tive systems to explore quantum coherence on a funda- 
mental basis and as possible basic unit {quantum-hit or 
qubit), for future quantum computers^ 2 -. 

The spin-boson model, where the TLS is bilinearly cou- 
pled to a bath of harmonic oscillators, is commonly used 
to quantitatively describe some aspects of the dissipa- 
tive dynamics of the above mentioned systems. In the 
spin-boson problem, it is common to evaluate the dynam- 
ics of the expectation value of the pseudo-spin operator 
(a z (t)) = P(t), as this is the quantity of interest in the 
experiments. 

To date, the spin-boson model has been mostly de- 
scribed within two main approximation schemes, each 
describing different regimes of temperature and coupling 
strength to the thermal bath. One main road of approx- 
imation is based on an expansion to leading order in the 
tunneling frequency A (see @ below), yielding the so- 
called noninteracting-blip approximation (NIBA^&i 3 -. 
Equations of motion for the dynamical quantity P(t) 
equivalent to the NIBA were also obtained using pro- 



jection operator technique s 14 ' 15 ' 16 ' 17 ' 18 . The NIBA equa- 
tions are easily solved numerically and yield analytical 
solutions in special cases^i 3 -. However, for a biased TLS 
they can be only used in the regime of high temperatures 
and/or strong friction. An improvement to the NIBA has 
been performed for a super-Ohmic bath, in the regime 
of strong coupling—. A more refined approach, which 
yields nonconvolutive dynamical equations, is the so- 
called interacting-blip chain approximation IBC A 20 ' 21 . 

In the range of weak TLS-bath coupling and low- 
temperature, however, the above mentioned schemes fail 
to properly describe the dynamics of a biased TLS, 
e.g. a symmetry breaking at zero temperature even 
for vanishing asymmetries is predicted. In this regime, 
path-integral method s 22 ! 23 ' 24 as well as perturbation 
theorie s 25 ' 26 are used. The lowest order perturbation the- 
ory and the path-integral approach have been shown to 
yield the same dynamics for weak Ohmic damping 2 ^. 

To date, only numerical technique s 28 ' 29 ' 30 ' 31 ' 32 can pro- 
vide a description of the TLS dynamics, being capa- 
ble to smoothly match between the high and the weak- 
coupling regimes. However, ab-initio calculations can be- 
come costly if the regime of low temperatures or the long 
time dynamics are investigated. 

In this work we propose a novel approximation 
scheme, which we call weakly-interacting blip approxi- 
mation (WIBA), capable to bridge between the weak- 
coupling and strong-coupling theories. It is based on 
the observation that bath-induced correlations among 
"blips" and between "blips" and "sojourns" are intrin- 
sically weak. Thus, in the WIBA, those correlations are 
included up to first order only. As in the NIBA, the time 
evolution of the population difference P(t) is given by 
the general master equation (GME) (t > 0) 



P(t) = - / dt'[K a (t-t')-W(t-t') + K s (t-t')P(t')] , 
Jo 

(1) 
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where the irreducible WIBA kernels K s and K a are, re- 
spectively, symmetric and antisymmetric with respect to 
an external bias e. As discussed below, such kernels are 
neither perturbative in the tunneling matrix A nor in the 
TLS-bath coupling, see Eqs. (J66J) , ([64]) and (|67j) . 

Since the inter-blip and blip-sojourn correlations be- 
come negligible at high temperatures or large system- 
bath coupling, the WIBA well matches the NIBA pre- 
dictions in this regime. On the other hand, at low tem- 
peratures and couplings such correlations, neglected in 
the NIBA, become essential to properly describe the dy- 
namics. In this limit, WIBA perfectly agrees with the 
predictions of the weak-coupling theories. Finally, in the 
intermediate coupling and temperature regime where the 
perturbative approaches fail, the WIBA yields a good 
agreement with predictions of ab-initio calculations. 

The paper is organized as follows: In Sec. |TT] we in- 
troduce the spin-boson model. Then, in Sec. IIIII the 
Feynman real-time path-integral formalism is shortly re- 
viewed, as it constitutes the starting point for the approx- 
imation schemes discussed in Sees. |IV] and [Vj Specifi- 
cally, we introduce the so-called Feynman- Vernon influ- 
ence functional and discuss an exact scries expression for 
the expectation value of the system "coordinate" P(t). 
Since the exact expression is nontrivial, the following 
Sections are dedicated to some approximation schemes. 
In Sec. IIV1 an extended version of the familiar non- 
interacting blip approximation (extended-NIBA) is pre- 
sented. There we derive a prescription to calculate the 
kernels entering an equation of the form (p}. Sec. [V] 
finally contains the major findings of our work. There 
we introduce the WIBA, which is able to bridge between 
the weak and strong coupling regimes within the same 
theory. The WIBA kernels entering ([1]) are found in 
analytical form. A comparison with predictions of ab- 
initio Quasi-Adiabatic Path-Integral (QUAPI) calcula- 
tions shows that the WIBA covers a wide spectrum of 
parameters and constitutes an interpolation between the 
NIBA and weak-coupling approximation schemes. 

II. THE SPIN-BOSON MODEL 

In this Section, the spin-boson model is shortly re- 
viewed, and the dynamical variables of interest are de- 
fined. To start, we consider the pseudo-spin Hamiltonian 
in the "localized basis" {\L), \R)}, 

H = |[^-Aa x ], (2) 

with a z and a x the Pauli matrices and hA the energy 
separation of the two levels at zero bias (e = 0) which 
accounts for the tunneling dynamics. Thus, we can inter- 
pret a z as a "position" operator such that (L\a z \L) = —1 
and (R\a z \R) = +1. 

We choose to model the environment as an ensemble 
of harmonic oscillator o 1 ' 2 ' 33 which linearly couple to the 
system "position coordinate" o~ z , ending up with a bath 



Hamiltonian Hb of the form (including also the interac- 
tion between bath and system) 

Hb = E 2 tar + m ^y] ~ E W'Xj • ( 3 ) 

3=1 3 3 = 1 

The whole system is thus described by the well known 
spin-boson Hamiltonian H = Ho + Hb- In the case of a 
thermal equilibrium bath, its influence on the system is 
fully characterized by the spectral density 

N C 2 

^=2^^-^ (4) 

3=1 3 3 

which reduces to a continuous spectral density once the 
number N of harmonic oscillators approaches infinity. 
Throughout this work, we choose spectral densities with 
power-law behavior at low frequencies, i.e. (s > 0) 

G{uj) = 2<W* 1 > s e-HM , ( 5 ) 

with 5 S being a dimensionless coupling parameter, o; p h 
a characteristic phonon frequency, and lo c the bath cut- 
off frequency which is taken to be the largest frequency 
in the model. Thus, this class encompasses the com- 
monly considered Ohmic spectrum with exponential cut- 
off G(uj) — 2au}e~\ UJ \/ UJc , with a = Si being the so called 
Kondo parameter of the TLS. 

In order to describe the system's dynamics, we focus on 
the time evolution of the reduced density matrix (RDM) 
p(t) = trBW(t), which we obtain after tracing out all 
bath degrees of freedom from the total density matrix 
W(i). We choose a factorized initial condition, namely 
at time t = the full density operator W(0) is expressed 
as a product of the initial system density operator (0) 
and the canonical bath density operator at temperature 
T. The initial preparation requires some attention^. One 
could distinguish two different preparations, according to 
the time when the coupling between system and bath is 
switched on, which we assume to happen at a time to < 0. 
In the first one, which we refer as "class A" , the bath is in 
canonical equilibrium and the system is being prepared 
at time to — in a certain state, e.g. corresponding to 
a z = +1. Then, the system evolves out of the state before 
the environment has relaxed to the shifted equilibrium 
distribution. This initial preparation is the typical situ- 
ation in electron transfer reactions when a specific elec- 
tronic donor is suddenly prepared by photoinjection2i. 
The other situation, "class B", is when the system has 
been held for a long time in a certain state, e.g. a z = +1, 
so that the bath had time enough to thermalize with the 
system. It corresponds to choose to — > — oo. Then, at 
time t = 0, the constraint is released and the system 
evolves with the spin-boson Hamiltonian. This initial 
preparation is performed, e.g. in rf-SQUID devices by a 
suitable choice of an external magnetic field. We shall 
refer ourselves to such initial preparation throughout the 
work. 
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The diagonal elements of the RDM are called popu- 
lations, whereas the off-diagonal terms coherences. In 
particular, we denote with 



sojourns 



PR® = Pa, 



l(t) = PRR{t) 



(6) 



the probability of finding the system in the right state at 
time t if it was prepared in the right state at time t = 
as well (with P^ one denotes the analogous quantity for 
the left well). 

In the spin-boson problem, it is common to evaluate 
the dynamics of the expectation value of the pseudo- 
spin operator (er 2 (i)) = pim(t) - phh{t) = P (t), namely 
the difference of populations in the localized basis. This 
quantity is also the quantity of interest of many experi- 
ments. Similar considerations as discussed in this work 
for P{i) can be done for the off-diagonal elements of the 
density matrix. 



III. REAL-TIME PATH-INTEGRAL 
APPROACH TO THE DYNAMICS 



In this Section we recall the main steps yielding an ex- 
act series expression for Pit), obtained within the path- 
integral approach. 



A. Influence functional and bath correlation 
function 

As we already anticipated, we assume a factorized ini- 
tial condition at time t — 0. For such initial condition, 
the exact formal solution for the RDM can be expressed 
in terms of a real time double path integral over piecewise 
constant forward cr(r) and backward <t'(t) spin paths^^ 
with values ±1. The effects of the environment are in an 
influence functional inducing non-local in time correla- 
tions between different path segments. Upon introduc- 
ing the linear combinations t](t) — [<t(t) + ct'(t)]/2, and 
£(t) = [ct(t) — £t'(t)]/2, the population difference reads 



P(t) = J D£D77.A[£,77]exp{$[£,77]}, 



(7) 



where A is the path weight in the absence of the bath 
coupling. A generic double path can now be visualized 
as a single path over the four-states of the reduced den- 
sity matrix, characterized by {ji(t) = ±f ,£(t) = 0) and 
(?7(t) = 0,£(t) = ±1). The time intervals spent in a 
diagonal (£(t) = 0) and off-diagonal (t](t) = 0) state are 
dubbed "sojourns" and "blips", respectively^ (see Fig. 
(H}). Due to the initial condition, the path sum runs over 
all paths with boundary conditions £(0) = £(t) = and 
?7(0) = I, n(t) = ±f . Environmental effects are included 
in the influence functional 



'0 



dhiih) S 2 ,i£(ti) + iR 2 ,ifi(ti) , (8) 



no = +1 









blips 



FIG. I: Generic path with 2n — 6 transitions at flip times 
ti, ti, ■ ■ ■ , t% n - The system is in an off-diagonal state (blip) 
of the reduced density matrix (RDM) in the time intervals 
t2j — (blue online) and in a diagonal state (sojourn) 

during the intervals t2j+i — t2j- Due to the initial preparation, 
the initial sojourn obeys the constraint 7/0 = +1. 



with the bath correlation function Q = S + iR being 

n , \ f 00 , G(uj) r , / fkj \. . . . 

Qit) = / auj — ^— cothl — — — (1 — coswi) + 1 smivt 
Jo L V2fc B J / 

(9) 

and Qj k : = Q{tj — tk)- For the class of spectral densities 
considered in Eq. ([5]), the functions S and R read^i 



S(t) = 2S s T(s - 1 



x cos [(s — I) arctan (u c t)) 




1 - (l + (w c t) 



2\ — 



2C(s- 1,1 + k) 



and 



R(t) = 2S s T(s 



1) f^Y 

\UpbJ 



(10) 



(11) 



x (l + (uj c t) 2 ) 2 sin [(s — I) arctan (uj c t)], 



where T(z) is the Euler's gamma function and £(q, z) 
is the Riemann's generalized zeta function^,. Moreover, 
k = l/hpLJ c becomes important as the ratio k^T/hw c 
becomes large. 

For a generic path with 2n transitions at times tj, 
j = 1,2, ...,2n, one finds £(t) = J2jZi^jH T ~~ *j) an( l 

f](r) = X)j=o 7 ?J^( T — H ere is Vo = 1 due to the 
initial preparation and £j = ±1, r]j = ±1 for j > 0. 
Because &j = —£,2j-i, the influence function in be- 

comes *W = $£L,b pa + *toL ( Fi S- ©■ The function 
describes intra-blip and blip-preceeding sojourn 



(n) 

intra. bps 



<F 

correlations, and reads 



intra, bps 



(12) 



$(") 4. $(") 

intra 1 bps 
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FIG. 2: Bath-induced non-local in time correlations among 
tunneling transitions. The interactions S2j,2j-i, and 
y,j_i (intra-dipole and blip-preceeding-sojourn interactions) 
which appear in the influence phase $i n tra,bps, cf. Eq. ((12}, 
are symbolized by the wiggled lines (blue and magenta online, 
respectively) . The double-dashed lines denote the inter-dipole 
interactions Aj t h, while the bold-dotted lines are the remain- 
ing blip-sojourn interactions Xj t k contained in the influence 
phase Winter, cf. Eq. ([16}. 



(n) 



n 



23,23-1 



^bps 1 



/] ZjVj-iYjj-i . 



where we split Xjj—i — i?2j,2j-i + Yjj-i, with 



Y 



3,3-1 



R2j-l,2j-2 — R 



2j,2j-2 



(13) 
(14) 

(15) 



(n) 

Moreover, the functional $ inter accounts for inter-blip 
and blip-sojourns interactions^ 



n j-l 

EE' 

0=2 fc=l 



= - E E ^Ai,k + * E E .v,. . (i6) 



n 3-2 

EE- 

j=2 fc=0 



The function A^/. contains the blip-blip interactions 
between the flip pairs {j, £;}, while the blip-sojourn in- 
teraction X; .fc yields a phase factor. To be definite, for 
k > 0, 

Aj,fc = S2j,2k-1 + S 2 j-X,2k ~ &2j,2k — S2j-l,2k-l , (17a) 
Xj,k= R2j,2k+1 + R%j-l,2k — R2j,2k — R2j-l,2k+l- (17b) 

The correlations -Xj,o depend on the initial prepara- 
tion, being of class A or B2. 



B. Exact series expression for P(t) 



with the conditional propagating function being 

rt / a2\» 



/•' / a 2 



X B n E ex P{ $( ™ } } 

{^=±1} {Vj=±l}' 



(19) 



with £?„ = cxp |— ie X)j=i £ j T 3 \> an d where we defined 

the blip length as tj = t^j — £2.3-1 • In (fl9|) . n counts the 
number of blips and the prime in {rjj = ±1}' means that 
the sum does not run over the initial and final sojourns, 
since they are fixed. Moreover, 



D{t 3 } = / dt 2n / dt 2n -i • • 



dti . (20) 



Performing the summation over the intermediate so- 
journs, one gets for the population differenced^: 



p (t) = 1 + E f^i} 

n=l Jo 



2 



(21) 



where B n s ^ 



= cos 



{eEU^3 



B, 



(a) 



,jTj | and 



n-l 

FW = G n JJ COS ((fo, n ), 
n-l 

= G n sin ((/>o,n) Y[ C0S i&k-n) , 



fe=l 



with <f> ktfl = J2]=k+i ^3 X 3,k and 
G n = exp {<He } 

{n n j—1 

- ^ S 2J,2j-l "EE && A i,* 
j=l j=2 fc=l 



(22a) 
(22b) 

(23) 
(24) 



The expression (j2"Tj) is still practically untractable. 
Thus, it is necessary to perform some approximations to 
describe the TLS dynamics. Two novel approximation 
schemes are discussed in the coming Sections IIVI and [V] 



The summation over the path histories reduces to an 
expansion in the number of tunneling transitions yield- 
ing a formally exact series expression for the population 
difference P(i)i*2. It reads (here we identify rj2 n with r)f) 

P(t)= E VfJ(vf,t;vo = i,o), (is) 

nj=±l 



IV. THE EXTENDED NON-INTERACTING 
BLIP APPROXIMATION 

In this Section we discuss an improvement to the 
familiar non-interacting blip approximation (NIBA)i>£, 
which better treats the blip-preceeding sojourns interac- 
tions. We call our more refined approximation scheme 
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A. Series expression within the extended-NIBA 



(b) 



t\ t2 t% ti ts t§ t 
Sztfj-l, R,2j,2j-1 



ti 



FIG. 3: Generic path and bath-induced correlations retained 
in the NIB A (a) and in the extended-NIBA (b). In both 
approximations the inter-blip and blip-sojourns correlations 
Kjjt and Xj t k : respectively, which appear in the influence 
phase Winter (see Eq. (JT3J), are neglected. Within the ex- 
tended-NIBA the blip-preceeding-sojourn interaction Yjj-i 
(magenta online), which contributes to the influence phase 
^bps (see Eq. (|14[0 and is being neglected in the NIBA, is 
retained. The first sojourn is treated differently, according 
to the initial preparation (here we choose the preparation of 
"class B"). 



" extended-NIBA" (see Fig. [3]). As we shall see, as the 
NIBA, the extended-NlBA enables to recast the series 
expression for P(t) into a generalized master equation 
(GME) of the form (JTJ with kernels of second order in 
the level splitting A. 

The advantage of the NIBA relies on its extreme sim- 
plicity and on the fact that it is non-perturbative in the 
coupling to the bath. Hence, the non-interacting blip 
approximation is a popular approximation scheme. Nev- 
ertheless, it has some intrinsic weaknesses, expecially in 
the asymmetric case (i.e. e ^ 0) for low temperature 
and weak coupling. For example, NIBA predicts the un- 
physical asymptotic limit ct^° n = — tanh(-^), implying a 
localization of the TLS (o"£° N = —1) at zero temperature 
even for vanishing asymmetries. 

The limits of validity of the theory are still dim, this 
approximation holding whenever the average time spent 
in a diagonal state (sojourn) (s) is very large compared 
to the average time spent in an off-diagonal state (blip) 
(t). Within the NIBA, the full inter-blip correlations 
Aj k and the blip-sojourn interactions Xj^ with j ^ fe+1 



(n) 



0). The blip-preceeding-sojourn 
To| are neglected as well (see 



{<Ke &£l a }. The explicit 



are neglected ($ 
interactions Yjj-i in Eq 

Fig. [3^). Hence, G n w exp 

form of the NIBA kernel is given in Appendix [A"l 

In the following Sec. [Vj we shall introduce a novel ap- 
proximation scheme, the weakly-interacting blip approx- 
imation (WIBA), capable to overcome these drawbacks. 
Before doing this, however, we need to introduce the ex- 



tended-NIBA where, as in NIBA, <fr 



(n) 



0, but the 



blip-preceeding-sojourn correlations i are retained, 
despite in approximate form. 



After performing the sum over the blip indices £j = ±1 



in Eq. (pH) , the extended-NIBA prescription <I> 
yields for the probability difference P e ^(t): 



(n) 



OO 

£ 

n=l 



P eN (t) = l + ^(-l)" / T>{tj}[g(Ti,so) 
Jo 







(25) 



+h(T 1 ,S )]'[[g(Tj,S j - 1 ) , 

where Tj = t%j — and Sj := £2.7+1 — t%j denotes the 

sojourn length. The extended-NIBA kernels are defined 
as 



g(rj,Sj-i) = A 2 e s(Tj) cos(eTj) cos[X,-j_i] 
h(Tj,8j-i) = A 2 e _5(T3) sin(er J )sin[A" : ,- J _i] . 



(26) 



Remember that = R(rj)— R(Tj +Sj-i)+R(sj-i). 

The functions <?(ti, sq) and /i(ri, sq) explicitly depend on 
the length sq of the initial sojourn, and assume a different 
form for preparation class A and B. 

In the following, we choose a factorized initial condi- 
tion at time t = with the particle being held at the site 
\R) [a z = +1) from time to = — oo till t = (class B), 
which amounts to consider so — > oo. Hence, Eq. (|26p for 
j = 1 reads 

gin, s ->• oo) = A 2 e" s(Tl) cos(eri) cos[i?(n)] = gn{n) , 

h(Ti,s -> oo) = A 2 e" s(Tl) sin(eri)sin[i?(ri)] = h N (n) , 

(27) 

being independent of sq- Notice that <?n(ti) and /in(t'i) 
coincide with the symmetric and antisymmetric NIBA 
kernels, respectively (cf. Eq. (|A1[) ). 

Due to the convolutive structure of Eq. (|2"5j) , it is eas- 
ier to evaluate the probability difference upon Laplace 
transformation. By exchanging the integration order and 
performing some change of variables, one gets 



1 r°° 
PeN(A) = - + £(-lf / -D 00 {r j ,s j :! 

n 

x [g(n,s ) + h(Ti,s )] Y[g(Tj,sj-i) , 

where P(A) = f Q °° dt exp [-Xt]P(t) = L t {P(t)} and 

poo poo poo poo 

/ Doo{Tj,Sj-i} = / ds n / dr n / ds n -i 
Jo Jo Jo Jo 

/•OO />oo />oo 

x... fox I d n ds e- x( -^7=^+ s ^+ s ^ 
Jo Jo Jo 



(28) 



ds n e 



(29) 



o 



I d£, n e~ Kn I dr n ... d£iC~ ACl / dn 
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with £j = Tj + Sj-i the length of each blip plus its pro- 
ceeding sojourn. Let us introduce the functions 

Ffo) = / d Tj g(r 3 ,^ - Tj ) , j > 1, (30a) 

JQ 

F fa)= ( 1 dn 5n(ti) , (30b) 
Jo 

Ao(Ci) = / l dn/»N(Ti). ( 30c ) 
Jo 

noticing that = i?(r,-) - Rfe) + - Tj). Then 

P(A) assumes the form 



oc n 

PcN(A)=-+-^n[-^ F fe)] w 

n=lj=2 

x{-% [*b(£i) + 4)«i)]}(A), 



(31) 



where the 1/A comes from the free last sojourn s n , or 
equivalently 



An(A) = - - -J2[-F(\)F-i[F (\) + A (\)} 

71—1 

1 lF (A)+i (A) 



A A l + F(X) 



(32) 



It is convenient to introduce the functions 
if e ° N (A) = Aio(A), i^ N (A) = AF(A) and 



W/ oN (A) = A F(A)-F (A) 



. Then Eq. 



becomes 



1 



P eN (A) = 



K^(X) - W / c n(A) 

A 

A + £* N (A) 



(33) 



B. Generalized master equation (GME) for the 
extended-NYQA model 



Eq. (|33p can be easily transformed back to the time 
domain. We hnd 



PMt) = - j dt'[K^(t-t') 

- W eN (t - t') + K s eN (t - t')P cN (t')} , 
with extended-NIBA kernels defined as 





d , . 


= Ht), 


(35a) 






:=A (t), 


(35b) 


Kits® = 


d , , 
dt F ^ 


■■=F (t), 


(35c) 


W*s(t) = 


Knit) - 


- K^ eN (t) . 


(35d) 



Although the time derivatives |35|) can be straightfor- 
wardly calculated from Eqs. (13"0"1) , the explicit depen- 
dence of the function g{Tj,^j — Tj) on £j still implies an 



(b) 



FIG. 4: Irreducible kernel K (s) (t) in the NIBA (a) and in the 
extended-NIBA (blue online) (b). The extended-NlBA kernel 
is obtained from the NIBA one by replacing the function R(t) 
with R(t) = R(t) —tR(t). Hence, the blip becomes "dressed". 



integral form for the function F(£j). An approximate 
form of F(£j) can be obtained if we observe that the 
average blip length (r) is suppressed by the intra-blip in- 
teraction exp [— S(t)] in the integrands of Eq. (|3TJ)) and 
that the imaginary part of the bath correlation function 
R(t), Eq. (fTTjl . slightly deviates from a constant. Hence, 
we approximate 



R(£j)+R(Zj-Tj 



-TjR(Zj) + 0[TfR(£j)}. (36) 



Corrections proportional to the second derivative of R(t) 
have been neglected. The functions <?(tj, Sj) and h(rj, Sj) 
defined in Eq. (j2"6"|) become (j > 1) 



9(Tj,Zj-Tj)ng(Tj,tj) 

:= A 2 e- S ^) cosfo) cos [R( Tj ) - Tj Rfo)] , 

H T i>Zi - t j)~ H T j^j) 

:= A 2 e- S ^) sin^r,) anfe) - TjRfa)] , 

(37) 

the corrections being of order 0[r?-R(^)]. In particular, 
we define <7eN(Ci) an d h e N(£j) as 

g*ifc) = = A 2 e- S fe) cos(s&) «»[£(&)] , 

haifa) = h(Zj,Zj) = A 2 c- 3 ^ sin( £ 0)sin[i?(0)] , 

(38) 

where R(t) = R(t) - tR(t). The extended-NIBA kernels, 
obtained as prescribed by Eqs. ([35]) and (|36|) . then read 



K, 



0,cN 



(*) 



= 9cN(t) , 
= hcN(t) , 

= Sn(*) , 

= 9eN(t) - gn{t) . 



(39a) 
(39b) 
(39c) 
(39d) 



The irreducible kernel K* N (t) entering the extended- 
NIBA master equation J33]) is shown in Fig. gb. The 
irreducible NIBA kernel (t) , Eq. (|Alh ) , is depicted in 
Fig. UK. As discussed in App. [A] within the NIBA an 
analogous GME as in (|34|) is obtained where, due to the 
approximation = 0, is Wn(£) = 0. 

The extended-NIBA kernels if| N and K° N differ from 
the conventional NIBA ones, and .KjJ, by the replac- 
ing of the imaginary part of the bath correlation function 
R(t) with the "dressed" one R(t) = R(t)-tR(t) (see Fig. 
HJd). At small coupling strength, the NIBA is recovered, 



since blip-preceeding-sojourn correlations become negli- 
gible. Despite its simplicity, however, the extended-NIBA 
already yields an improvement to the NIBA in the inter- 
mediate coupling regime (see Fig. [5]). 
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FIG. 5: A comparison between the standard NIBA and the 
extended-NIBA for an Ohmic bath (s = 1) is shown. The 
parameters are: T = 0.1, a = 0.1, e = 1 (in units of the 
tunneling frequency A). One can see that in the chosen inter- 
mediate parameter regime the extended-NIBA coincides with 
the conventional one at short times. However, it predicts a 
different asymptote from the NIBA one (see inset). Hence, 
at moderate damping and temperatures, the blip-preceeding- 
sojourn correlations retained in the extended-'NIBA become 
important. 



V. THE WEAKLY-INTERACTING BLIP 
APPROXIMATION (WIBA) 

To bridge between the strong damping situation de- 
scribed by the NIBA and the extremely underdamped 
case we observe that, for spectral densities of the form 
(O, the blip-blip interaction terms Aj t k as well as the 
blip-sojourn terms Xj t k (k j '< — 1) are intrinsically 
small. Therefore, we propose a novel approximation 
scheme, which we call weakly-interacting blip approxi- 
mation (WIBA). Within the WIBA, the full $1™^^ is 
retained as in the extended-NIBA and one expands the 
influence functional exp {$j" ter } up to linear order in the 
blip-blip and blip-precccding sojourns interactions Aj } k 
and Xj,k (see Fig. [6|). Hence, 



exp{$ (n) } m exp{$[ 



(«) 

intra, bps 



>( 



l + $ . 



(40) 



In other terms, all contributions which involve 4>k,n (cf- 
Eq. (|22|) ) must be expanded up to first order in Xj t k, j > 



2, and G n « exp <^ 5He 



(n) 



}( 



(») 



" intra, bps 

(cf. Eq. [[23])), All terms of the order AX or higher 
are neglected. As usual, the first sojourn and blip must 
be treated differently from the others, according to the 
initial preparation. 



x v S2j,2j-1,-R.2i,2j-1 




FIG. 6: Generic paths and bath-induced correlations retained 
in the WIBA. The linearized influence functional <&[™ ter yields 
one single blip-blip correlation (a) or blip-sojourns interaction 
(b) for each path. 



The inter-blip and blip-sojourns correlations Aj t k and 
Xjjf are known to become essential to properly describe 
the dynamics of a biased TLS at low temperature a 22 i 23 ' 24 . 
Specifically, in Ref. [23| a systematic weak-coupling ap- 
proximation (WCA) was developed where all bath- 
induced correlations were linearized: exp 

It was shown in Ref. H3 that the WCA exactly 
matches results obtained within the lowest-order Born 
approximation for weak-coupling to an Ohmic bath. 

Hence, the WIBA constitutes an improvement of the 
extended-NIBA model on one side and of the WCA on 
the other. At high temperatures, where the dipole-dipole 
interactions are negligible, the WIBA kernels reduce to 
the extended-NIBA ones. By expanding the WIBA ker- 
nels to first order in 8 S , the weak damping kernels in Refs. 
22.23j are recovered. 

In contrast to the extended-NB3A or the weak-coupling 
approximation, the WIBA has no small parameter in a 
strict sense, and it is based on the consideration that the 
correlations $j" ter are intrinsically weak over the whole 
parameter regime. As we show below, the WIBA indeed 
well describes the TLS dynamics over a wide parameter 
range. 



A. Series expression within the WIBA 

Let us again start from Eq. (|2"Tj) and apply the WIBA 
prescription. After performing the sum over £j = ±1 we 
then find 

PwibaW = 1 + / dt 2 f dti[-{ gi + hi)} 
Jo Jo 

co „t ( In 

- I 1 ' 1 \ 



n I n 



1=2 \ j=l+l 



l-l 



h \ X\gj\ {gi + hi)X lfi 



8 



n I— I 

•EE I n* 

;=3 fc=2 \j'=i+l 



/-i 

+ hi | J~[ | hkXi ik -i 

\3=k+l 
=P eN (t)+ Pinter (*) 



ft/ 1 n ft i ftfcA/,fc 
fc-i 

n » i Cffi + fto 



(41) 
(42) 



where P e N (t) is the exienrfed-NIB A expression (|25|) , while 
Pnter (t) contains the contributions coming from the lin- 
earized influence functional exp{$^J? er }, cf. (|3D|). A 
graphical representation of the time-ordered sequences 
entering Eq. (|H|) is shown in Fig. [6l In (|4l"jl . <?j and ftj 
are the functions already occurring in 



(43) 



fJj =9(Tj,Sj-l) , 

hj = ft(r,-,S3_i) . 

The functions g., and ftj are analogously defined as 

9i = AV 5(t ^ cos(erj) sin^-j-x] , 
hj = A 2 e~ 5(Tj) sin(£Tj-) cos[JCjj-_i] . 



(44) 



Notice that gj and ^ are symmetric in the bias, while hj 
and hj are antisymmetric. Moreover, 



9i = 0n(ti) , 
fti = ftN(n) , 



(45) 
(46) 



cf. Eq. ([27]). and analogously for g\ and hi. To proceed, 
let us consider the function (for k > 1) 



£i ( - fe+1 >(i 2/ -i 2fe - 2 ) = 



dt 2 i-i 



1 2k 



dt 



2k-l 



l-l 



(47) 



< ft/ I Yl 3j ] [ftfcAf.fc + ftfe-X/,k-i] 



entering (JUJ), where 2(2 — fc + 1) denotes the number 
of tunneling transitions. Upon introducing the variables 
£ = t%i — t 2 fe-2 and £j = 7j + Sj—i, Eq. (|47p assumes the 
form 

£^ +1) (0 = jf d$ k J k d£ lThl (&)/(£ - 6 - eo (48) 

x [7h* +7h„(6)A/,fc] , 

where 7^ . ^ . (£j ) is an operator which acts on a generic 
function w(jj) as (e.g. let us consider the operator 

lh 3 (&)) 

7fej (0) M := / drj /i (73, & - Tj)w{Tj) . (49) 
Jo 

In (148|) . the function / is given by the expression 
/(*2/-2 — ^2fc) = 8{t 2 i- 2 — hk) 



^21-2 r^2k + 2 

+ / dt 2 l-3 ... 



l-l 



dt 2k +i [ n ~9j ] • (50) 

i=fe+l 



Upon introducing s = t 2 i-2 — t 2 k and £j = Tj + Sj-\ 
as in Eq. ESl. we find 



/(«) = 8(s) 



dik+i 

n JO 

i-3 



£fc4 



rfTfc+i [-g(r fc+ i,£ fc+ i - Tfc+i)] 



^6-2 / rfTi-2 [-5(71_2, 6-2 - T;_ 2 )] 

Jo 



dri-i 




(51) 



where the first derivative of the conditional probability 
jJ e N has been introduced. It satisfies the master equation 
for the extended-NIBA with the symmetric kernel only: 



PeN(t) = - / dt'K s eN (t - Op-n(0 ■ ( 52 ) 
Jo 



Analogously, we can treat the contributions to Pi n t cr 
which depend on the initial preparation. Specifically, we 
introduce for the case k = 1 (here £ = t 2 i — 0) 



£8,(0 = I «i / dei-m, (&)/(£ - a - £0 

Jo Jo 

x [731(6)^,0-79! (6 )A/,i] , 

£8,(0 = / ^ / l <%ir hl mm - ti - 6) 

Jo Jo 

x [7fci(fi)-Xi,o+Tft 1 (6)Aj,i] , 



s 



(53) 



(54) 



where the operators 7 Sj /g . (£,■ ) act analogously as in (l4*9"j) , 
and/( S ) = «5( S )+^ N 2) ( S ). 

As discussed in the Sec. IIV B[ we choose a factorized 
initial condition at time t = 0, which corresponds to set 
so — ► 00. This means that in Eqs. (|53p and (|54[) . A^o = 

P(*2Z ~ *l) ~ P(*2/-l - *l) = P(£ " £l + ri) - P(£ - £l + 

n - 77). Moreover, A z .i = S(f - £1 + n) + S*(£ - £1 - 

7l)-S(£-£l)-S«-£l+71-7)). 

As before, it is more convenient to work now in the 
Laplace space, where many terms factorize. In fact, 
one can identify (see e.g. Fig. [5]) products of one or 
several irreducible extended-NIBA kernels (P(A), A (X) 
and Fq(\)) with one of the irreducible kernels contain- 
ing inter-blip and blip-sojourns interactions (Sg (A), 
Eq2(A) and Sqs(A)). After exchanging the integration 



order as done in Eq. 



the probability difference reads 



9 



in the Laplace space 



Pwiba(A) = P eN (A) 



1 



GO I 71 

EE 

' n=2 I 1=2 
n l—l 



n - p w 

3=1+1 



EE 



Z=3 fc=2 
fc-1 



n 



n 



^0 : a + Sq ;S 



(J-k+l) 



(55) 



x 



n - p w 

3=2 



(f (X)+A (X) 



After some changes of variables, Eq. (|55|) can be recast 
in the form 



Pwiba(A) =P eN (A) 



(7+2) 



-F(A) 



1 oo m ( 

+1EE [S0.. + E 
m=0 7=0 I 

r „ 1 (7+2) r » -i m-7 

+ [S S J -F(A)J (m-7 + 1) 
(F (A)+i (A) 



m— 7 



(56) 



which becomes, after noticing that X!m=o 

7—0 Z_-/m=7 ' 



Aviba(A) = F cN (A) 

+ TE [^a + So,,] TRW 
m=0 I ~ p=0 



(m+2) 



oo 

]T[-A(A)] P ( P + 1 



p=0 



F (A) + i (A) 



(57) 



It is convenient to introduce the irreducible kernels 



E » = Em=o E * and S o,a/ s = E„=oV ob- 
tained upon summing over the number of tunneling tran- 
sitions. Thus, recalling the extended-NIBA expression 
([281). we find 



^oo A(m+2) 



Aviba(A) 



1 1 

A + A 



(i?o(A) + i (A) 
x ]T [-F(A)]" 1_1 [-F(A) + mS s (A) 

77) =0 

+ - [So, a (A)+E Q , s (A)] £ [-F(A) 



(58) 



Keeping in mind that we should always retain the in- 
teractions A and X up to the first order, Eq. can be 



rearranged in a more compact form as 
1 1 
A + A 



Pi 



WIBA 



(A) 



F (X)+A (X) 



do 

E , a (A) + S ,s(A)] £ [-F(A) + E,(A) 



m=0 



(59) 



yielding 
Aviba(A) 



1 



(f (A) + i (A)) + E , o (A) + So,. (A) 
1 + F(A) -E S (A) 
1 + Wwiba/A - F W iba/A 



a + Ar^, IBA 

where we introduced the functions 



A, 



(A)=A F(A)-E g (A) 



WIBAV 

^wiba(A) = A [i (A) - S , a (A) 
^oV I ba(A)=a[ j F (A)-S 0;S (A)] , 
Wwiba(A) = A^ IBA (A) - Ap jWIBA (A) . 



(60) 
(61) 

(62a) 
(62b) 
(62c) 
(62d) 



B. Generalized master equation (GME) for the 
WIBA model 



After multiplying both sides of Eq. (|61| by A+ A| VIBA , 
one can recognize the Laplace transform of the general- 
ized master equation Eq. (TT]), where from Eq. (|62[) the 
WIBA kernels in the time domain are 



A, 



WIBA 



A, 



(t) = F(t)-E.(t), 
{t) = F (t)-t 0tS (t), 
K^ 1BA (t)=A (t)-± , a (t), 



0,WIBA 



(63) 



WIBA 



(*) — ^WIBAW ~ ^0, WIBAV*) • 



The prescription (|63[) allows the explicit evaluation of 
the irreducible kernels. This procedure is illustrated in 
App. El Explicitly, the antisymmetric WIBA kernel 
reads 

#wiba(0 - Kn(0 

S /•£-& 

o Jo 



x \m(ti) x f,o -5n(6)A/, 



1J • 



(64) 



where, as in the extended-NIBA case, cf. (j3"6")l . correc- 
tions of order have been neglected. The func- 
tions <?eNj ft-eN have been defined in {38]). Analogously, 
we also introduced here the functions 



3cn(£j) 



A 2 e~ s ^ cos(eQsm[R(Q] , 
A 2 e - s fe'sin(£^)cos[^(^)] , 



(65) 
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FIG. 7: Irreducible symmetric WIBA kernel. Notice that the 
lowest order reproduces the extended-~NTBA case, see Fig. fib. 
The inner bubble represents the contribution coming from the 
sum of all orders in A 2 which coincides with the symmetric 
part of the solution of the master equation for P(t) in the 
extended-NIBA case. 

with R(t) = R(t) - tR(t). The symmetric WIBA kernel 
becomes 

*W(0 = ^e N (0 

d£l C^//leN(C/)PeN(C-^l ~ £/) 

Jo Jo 

x [ft. cN (6) A /,i] ■ 

(66) 

The irreducible symmetric kernel is shown in Fig. [7] 
Moreover, the function IFwiba is now given by 

Wwiba(0 = ^wiba(0 - #o,wiba(£) , (67) 

where 

^0,WIBa(£) = ^0,eN(S) 

- / d^i / d^ f h eN ^ f ) PeN ^-^ -£ f ) 
Jo Jo 

x [ft N (a)x / , +ft N (a)A / , 1 ] . 

(68) 

In the Eqs. and (|B8|), X/, = R(£) - R(£ - £/), 

and for all the three kernels A/ 1 = — 5(£ — £i) + 

5(^-^-0)-^-^/). 

C. Dynamics within the WIBA 

Let us now come to the WIBA predictions. In this Sec- 
tion we will be showing the comparison of all the theories 
previously discussed, namely the NIBA, its extended ver- 
sion, the weak-coupling approximation (WCA) and our 
WIBA, for different choices of the parameters. In the fol- 
lowing subsections we examine the behavior of P(t) in the 
case of Ohmic bath (s = 1 in the spectral density |(SJ)), 
with cutoff frequency to c — 50A, and of super-Ohmic 
dissipative environment (we choose s = 3), for different 
choices of the frequencies oj c and w p h- 

1. Ohmic case (s — 1) 

We expect WIBA to work particularly well in the 
Ohmic case due to S(t) ~ t at long times. This has the 



simultaneous effect of suppressing the blip lengths and to 
yield a vanishing interblip interaction at large blip sepa- 
ration. 
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FIG. 8: Time evolution of P(t) for an Ohmic symmetric two- 
state system. The WIBA (full lines) coincides with the NIBA 
(dashed lines) and the extended-'NIBA (dotted lines) over the 
whole range of parameters. Here a = 0.1 and T = 0.1 (in 
units of A) have been chosen. For such a choice, the WCA 
(dot-dashed lines) slightly deviates from all other predictions. 

We first present (Fig. [5]) the results for an Ohmic sym- 
metric two-state system (e = 0) for a generic choice of 
the parameters. For the symmetric case, the NIBA is 
expected to predict the correct time evolution over the 
whole range of coupling strengths and temperatures. One 
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FIG. 9: Time evolution of P(t) for Ohmic damping and finite 
bias. The chosen parameters are a = 0.01, T — 0.1 and 
e = 1 (in units of A). The WIBA (full lines) coincides with 
the WCA (dot-dashed lines) for low temperatures and small 
coupling, while the NIBA (dashed lines) and the extended- 
NIBA (dotted lines) predict an unphysical asymptote, for an 
asymmetric system. 

sees a complete agreement of the WIBA with the {ex- 
tended-) NIBA. For the chosen parameters, also WCA 
well agrees with the WIBA predictions. 

The situation becomes more intricate for the case of 
finite bias (e ^ 0), since on the one hand NIBA is ex- 
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pected to fail at low temperatures and damping, while 
the WCA becomes inappropriate at large temperatures 
and/or damping. In the following, we fix the external 
bias as e = A. 

Fig. [9] shows a case of low temperature k^T < E = 
\/e 2 + A 2 and small damping a, namely a — 0.01 and 
T = 0.1 (in units of A). One can see that the agreement 
between WIBA and WCA is striking, whereas both NIBA 
and extended-NIBA fail to reach the correct asymptotic 
value, predicting an unphysical symmetry breaking. 
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FIG. 10: Time evolution of P(t) for a = 0.1, T = 0.1. The 
WIBA (full lines) exhibits more coherence than the QUAPI 
(lines with bullets), being however its predictions closer to 
the QUAPI than the WCA (dot-dashed lines). The extended- 
NIBA (dotted lines) is also getting closer to the QUAPI pre- 
dictions, the NIBA (dashed lines) still predicting a strong 
localization in the left well. 

In Figs. [TU1 HU and a comparison between 

WIBA and the numerical ab-initio path-integral ap- 
proach (QUAPI ? ?) is made, in order to proof the validity 
of our theory in a regime of intermediate-to-high temper- 
atures and coupling strength. As the coupling strength 
is raised to a = 0.1 (Fig. [Tuj) . while keeping the tem- 
perature constant, one sees that the asymptote of the 
extended-NIBA is slowly moving from the NIBA one to- 
wards the QUAPI one. At the same time, the WCA 
also disagrees with the WIBA predictions, the former 
starting to be invalid for higher coupling strength. The 
WIBA, however, shows the presence of some unphysi- 
cal "beatings" at the onset of the long-time dynamics, 
whose origin is not yet well understood. Hence, some 
more investigations seem to be required, in order to bet- 
ter understand the theory in this regime of parameters. 

As one raises the coupling strength (a = 0.25) (Fig. 
ITT1) . by lowering the temperature (T = 0.01), one can 
observe that the WCA completely fails in predicting 
the right asymptote or the intermediate dynamics. The 
NIBA also predicts a wrong asymptotic value, while cor- 
rectly describing the short-time dynamics. For such a 
choice of parameters, the extended-NIBA very well agrees 
with QUAPI, although some oscillations are still present 
in the theory. The WIBA also smoothly oscillates close 
to the QUAPI predictions. 
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FIG. 11: Time evolution of P(t) for a = 0.25, T = 0.01. 
The WCA (dot-dashed lines) completely fails in describing 
the short-time dynamics and also predicts a wrong asymptote. 
The WIBA (full lines) as well as the extended-NIBA (dotted 
lines) very well agree with the QUAPI predictions (lines with 
bullets), despite some spurious oscillations, still contained in 
the theory. The disagreement with the NIBA is striking, al- 
ready at this intermediate-to-high regime of parameters, al- 
though the NIBA correctly describes the short-time dynamics. 



Finally, Fig. [12] shows that, for the parameters a = 
0.25, T = 0.1, the WIBA and the extended-NIBA are al- 
most indistinguishable from the QUAPI predictions. The 
WCA and the NIBA fail to describe the dynamics, even 
though the NIBA correctly reproduces the short-time dy- 
namics. By looking at Figs. [TUl[TT1 and[T2l it is interesting 
to notice that the WIBA predicts more coherence than 
QUAPI. The overall agreement with QUAPI though re- 
mains very good. 

Hence, we showed that WIBA works very well for 
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FIG. 12: Time evolution of P(t) for a = 0.25, T = 0.1. In 
this regime the extended-NIBA (dotted lines) as well as the 
WIBA (full lines) are almost indistinguishable from QUAPI 
(lines with bullets), whereas both WCA (dot-dashed lines) 
and NIBA (dashed lines) fail in describing the dynamics, the 
latter being, however, able to reproduce the short-time dy- 
namics. 
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strong damping and/or high temperatures and small cou- 
pling strength and/or low temperatures, where is able to 
reproduce previous perturbative theories. In an interme- 
diate regime of parameters, WIBA is shown to give an 
appropriate description of the dynamics. 



On the other hand, Fig. [14] shows the WIBA predic- 
tions for P(t) when the cutoff frequency lo c is of the or- 
der of the tunneling frequency A, namely when the bath 
becomes "slow", a case where the WCA and the NIBA 
are expected to fail. This case is the most difficult one, 



2. Super- Ohmic case (s = 3) 

Let us now consider the predictions of the WIBA in 
presence of a super-Ohmic bath. We choose to show here 
a very common situation, corresponding to s = 3. In con- 
trast to the Ohmic case, analyzed above, for the super- 
Ohmic case is unclear if the WIBA can perform well over 
a broad regime of parameters. In fact, since S(t ~ A -1 ) 
differs only little from its asymptotic value S(t — > oo), 
the interblip interactions are weak and the WCA is ex- 
pected to be a good approximation in a wide regime of 
parameters. This also implies that S(t) is not effective 
in suppressing long-blip lenghts and the NIBA might not 
be justified. Notice that we expect WIBA to work well 
in all the situations such that 1 < s < 2, since for those 
parameters the function S(t) increases with time. 

In the following figures, we keep temperature, exter- 
nal bias and coupling constant fixed (T = 0.1, e = 1, 
5% = 0.01, in units of A), while keeping the ratio Lo c /u) p h 
constant (o; c /a;ph = 8). This choice is reasonable and 
agrees with some experiments^. The ratio w c /w p h acts 
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FIG. 13: Time evolution of P(t) for 8 a = 0.01, u c = 200 
(T — 0.1, e = 1). Here one sees no difference between WCA 
(dot-dashed lines) and WIBA (full lines), showing that the 
chosen parameters lie in an effective weak-coupling regime. 

as an effective correction factor to the real coupling 
strength, see Eqs. (jTUJ) and (fTTj) . Moreover, since we 
fix the temperature, in the bath correlation function cor- 
rections proportional to the temperature (k = l/h(3uj c ) 
become important as the cutoff frequency is decreased. 

In Fig. [13] the cutoff frequency uj c = 200 is assumed to 
be very large compared to the system frequency scales. 
There, one sees that the WIBA well matches the WCA 
predictions, as expected in the regime of weak damping. 
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FIG. 14: Time evolution of P(t) for S s = 0.01, u) c = 1 
(T = 0.1, e = 1). Here we notice discrepancies between WIBA 
(full lines) and QUAPI (full lines with bullets), since the ex- 
tended-NTBA predictions (dotted lines) also account for too 
many oscillations. Nevertheless, it gives much better results 
than NIBA (dashed lines), predicting unphysical values for 
Pit). Finally, the WCA (dot-dashed lines) is so far the best 
model in this regime, although it lies apart as well from the 
QUAPI predictions, since the perturbative approach for such 
parameters begins to fail. 



since the bath is very coherent and memory effects are to 
be taken into account, which requires to perform a very 
good description of the full bath dynamics. One sees 
that the NIBA completely fails to reproduce the dynam- 
ics, even reaching unphysical values. The extended-NIBA 
works better, approaching closer the QUAPI predictions. 
Nevertheless, too few correlations are taken into account, 
and it oscillates still too much with respect to the nu- 
merical plot of QUAPI. The WIBA shows discrepancies 
from the QUAPI as well, being still "too" coherent, even 
though its predictions are more accurate than the ex- 
tended-NlBA. The WCA, despite better than WIBA in 
this regime, also lie apart from the numerical prediction 
of QUAPt. In this range of parameters, the temperature- 
dependent corrections in the bath correlation function 
Q become relevant and the perturbative weak-coupling 
approach begins to fail. Thus, more analysis of the com- 
plicated super-Ohmic case is to be done, in order to bet- 
ter understand the different dynamical situations which 
take place by varying the coupling strength 5 S , the cutoff 
frequency lo c and the phonon frequency a; p h- 
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VI. CONCLUSIONS 

To conclude, we presented here a novel analytical 
scheme dubbed WIBA (Weakly-Interacting Blip Approx- 
imation) which is able to match between diverse approx- 
imation schemes for different parameters choices. In 
particular, the WIBA is valid over a large regime of 
temperatures and coupling strength for a wide class of 
spectral densities ([5]). The WIBA reproduces very well 
weak-coupling approximation schemes for small damping 
and temperature and the well-known NIBA in the oppo- 
site regime of strong damping and high temperature. It 
yields a good, though not perfect agreement, with ab- 
initio QUAPI calculation in the regime of intermediate 
temperatures and damping. Hence, the WIBA overcomes 
the dim validity of the perturbative approaches discussed 
in the previous Sections which, up to now, was only pos- 
sible via numerical ab-initio schemes. It also serves as 
benchmark to limiting simple analytical schemes, in or- 
der to proof their range of applicability over the various 
regimes, and combine weak-coupling and strong-damping 
approaches in one single, unified model. 

Our approach is based on the consideration that bath- 
induced blip-blip and blip-sojourn interactions are weak 
in the whole regime of parameters for the spectral densi- 
ties ([5]). It is expecially needed in several contexts, and 
we mention here few of them: i) The WIBA could be 
very useful when experimentally the bath temperature 
or the external bias are varied over a wide range: In fact, 
WCA and NIBA are unreliable at high temperatures and 
low-to-intermediate driving, respectively, ii) It should be 
used to investigate the situation of several TLS's inter- 
acting with a common heat bath^l, as in glasses^, at low 
temperature being characterized by a broad distribution 
of tunneling parameters and asymmetries. 

We must, however, notice, that in the case of "slow" 
environments with cut-off frequency of the order of the 
tunneling frequency, still our model has to be improved, 
since neither the WIBA nor other analytical approxima- 
tion schemes are able to reproduce the correct onset of 
decoherence which in fact takes place. This situation oc- 
curs e.g. in non-adiabatic electron transfer— or for charge 
qubits interacting with piezo-electric phonons^. 
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APPENDIX A: DERIVATION OF THE NIBA 
KERNELS 

In the non-interacting blip approximation, as men- 
tioned in the Sec. IIV[ one neglects all blip-blip cor- 
relations and blip-sojourns interactions in Eq. (|16p . 
i.e. 3>i"ter = 0- The blip-preceeding-sojourn interac- 
tions are neglected as well (3>j5,i>s = 0), i.e. one sets 
Xjj-i s» i?2j,2j-i- As a consequence, the influence func- 
tion (fT2"]) factorizes into individual influence factors de- 
pending only on the dipole length Tj := t%j — i2j_i- Fol- 
lowing an analogous procedure as in Sec. IIV A[ the NIBA 
series expression can be derived, and the corresponding 
GME, of the same form as in Eq. is obtained. The 

kernels finally read 

K^{t) = A 2 e~ S(t) cos(rt) cos[R{t)} , (Ala) 

K$(t) = A 2 e- S(t) sin(et) wa[R(t)] , (Alb) 

Kfa(t) = A 2 e~ s(t) cos(ei) cas[R(t)] , (Ale) 

Ws{t) = K^t) - Kfa(t) = . (Aid) 

The kernels are of lowest order in the tunneling matrix 
A but are non-perturbative in 5 S . The NIBA is expected 
to be a good approximation whenever the average time 
spent in a blip is much larger than the time spent in a 
sojourn. 

In general, the inter-dipole interaction can be safely 
neglected for all temperatures for sub-Ohmic damping 
s < 1, and the TLS is expected to exhibit incoherent 
dynamics even for very small coupling 8 S . For Ohmic 
and super-Ohmic damping, NIBA is expected to be a 
good approximation at high enough temperature and/or 
strong damping. In particular, spectral densities of 
the Ohmic form reach for large temperatures faster the 
asymptotic behavior S(t) oc t at long times, implying 
that the correlations in Aj^ cancel out exactly. 

APPENDIX B: EVALUATION OF THE WIBA 
KERNELS 

In this Appendix, the prescription ([63]) to evaluate the 
WIBA kernels in the time domain is illustrated. Let us 
consider, to fix the ideas, S , a (0 = Y^i a(Q- The 
kernel Eq ' , given in Eq. (|53p . is reported here for clarity: 

= J #i / d&TS, - & - 6) (B1) 

x [7 Sl (6)A;, -7si(Ci)A;,i] , 

with X lfi = Rfc - a + n) - R{£, - & + n - n). Notice 
that the operators 7's depend on the functions g/g (h/h) 
defined in Eqs. flU} and (05]). 

In the function f(s) which appears in the expression for 
the kernel (see Eq. ([50]) ). the first derivative of the con- 

(1—2) 

ditional probability p cN (s) is present. Then, it could 
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be integrated out by parts, obtaining 

ri-£i 



x P C ' N % - Si ~ SO [7« (Si)X,o - 7 fll (Si)A U ] • 

(B2) 

Here, the first derivative of the operator 7/j ( (Sz) acts on 
the functions Xi_ = i?(S - Si + n) - -R(S - Si + T i - tj) 

and A;,! = s(f - a + n) + S(£ - 6 - n) - S(£ - 6) - 
S(S - a + n - n ) as 

7&(6)*i,o 

[tl _ 

= % / dn h(n,£i - Ti)Xifi 

(B3) 



% / dTih(Ti,£i)Xi t0 
Jo 

MSO [R(t - Si + r x ) - - 6 + n - 6)] , 



and 



7h(6)A U 
/■& _ 

= % / dri h(n,^i - n)A lA 



J o 

= WSO [S(£ - 6 + n) + S(£ - & - CO 

As already observed when calculating the extended- 
NIBA kernels, the explicit dependence of the integrand 
of Eq. (|B3p on yields an integral form for 7x(S0^z,o- 
An approximate form is obtained if we apply the same 
prescription (|36[) as for the extended-NIBA. Namely, one 
approximates the bath correlation differences in Yjj—i 
with the first derivative of R(t). Hence, for g and h Eq. 
(|37| holds, whereas for g and h it follows (/ > 1) 



g(n,& - r;) w #(t;,S0 

= A 2 e~ 5 ( Ti ) cos( £ r i ) sin[i?(r i ) - r^(6)] , 

= A 2 e" s ( Ti ) sin( £ rO cosfiZfa) - 7j.R(£i)] , 

(B5) 

the corrections being of order 0[rf R{^i)\. We define 

<7eNte) = §($i,£i) and /i oN (S0 = ^(ShSO- 

In order to evaluate the antisymmetric WIBA kernel 
-^wiba fr° m the prescription (|63[) . one sees that the pre- 
vious expression must be derived with respect to £. After 
substituting u = £ — Si in the Eq. {B2| , we get 

4J(S) = / <*« fW^e N (S0 

Jo Jo 

X PcN 2) ( w ~ SO frsi (S ~ «)-Xi,o - 7si (S - «)Ai,i] , 

(B6) 



with X; o = R(u + ri) - i?(w +rx - SO an( i A; i = S(u + 
n) + S(u - &) - S(u) - S(u + n - SO- 

Then, in order to get the derivative of £q „(S) with 
respect to £, it is enough to evaluate 

hgx (S - «)^,0 ~ 7si (S ~ w)Ai,i] 
/•£-« 

= c\ / dn [5i(n,s -> oo)A;. - gi(rx,8 -» oo)A ; .i] 
Jo 

= 5n(S - fi)A; !0 - 5n(S - w)A;,i , 

(B7) 

where now X t = R(£) - - £0 and A; i = S(Q + 

S(S-Si-SO-S(S-Si)-S(S-SO- 

Hence, we obtain for the Z-th order (after the deriva- 
tive, we substitute back £i = £ — u) 

^£(S) 

= f <%x t d& MSOpIn 2) (S - 6 - CO ( B8 ) 
Jo Jo 

x [5n(6)A;, -5n(6)A(,i] • 
In the end, the antisymmetric WIBA kernel reads 



(B4) = ^cn(0 - £ A^ IBA ( ff) (C) 



-/ <%1 / ^//ieN(C/)PeN(C-^l -C/) 
Jo JO 

x [<7n(Si)^/,o -5n(6)A/,i] , 



(B9) 



where we introduced PcN (^) = X)^L2PeN ^(O- Following 
similar lines, we obtain for the symmetric WIBA kernels: 

oo 

K S , W1BA (0 ee Kl cN (0 - J2 ^o,wiba (<t) «) 

(BIO) 



(T=2 



d£l / ^/ /leN (0 ) PeN (C ~ Si ~ ) 



JO 



and 



x [ft N (Ci)jr/, + /i N (a)A/,i] , 

oc 

= ^ C S N (0 



cr=2 



i rt-Si 
d£x / C?S/ h eN (S/ ) PeN (S - Si - ) 



JO 



(Bll) 



x [ft e N(Sl)A/,: 



We identified in all equations the last time interval with 
Tf and correspondingly £/. In this last equation, we also 
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identified the first time interval Tk with n (therefore £fc 
with £1) and Af & with A/ i. 

In Eqs. JB9| and fBTO| '. is X ffi = R(£) - R(£ - £/), 
and for all the three kernels A/,i = fi^ + S^— £1— £/) — 
S{€-€i)-S(£-£f)- InEq. dlll|), the term proportional 



to Xf.o is absent because it vanishes after the integration 
by parts and the derivative, since it does not depend on 
s . 
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